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ABSTRACT 



Average charged multiplicities have been measured separately in b, c and light quark 
(u, d, s) events from Z° decays measured in the SLD experiment. Impact param- 
eters of charged tracks were used to select enriched samples of b and light quark 
events, and reconstructed charmed mesons were used to select c quark events. We 
measured the charged multiplicities: n uds = 20.21 ± 0.10 (stat.) ± 0.22 (syst.), n c = 
21.28 ± 0.46 (stat.)laM (syst.) and n b = 23.14 ± 0.10 (stat.)!":!? (syst.), from which 
we derived the differences between the total average charged multiplicities of c or 
b quark events and light quark events: An c = 1.07 ± 0.47 (stat.)!^ (syst.) and 
An b = 2.93 ± 0.14 (stat.)!^ (syst.). We compared these measurements with those 
at lower center-of-mass energies and with perturbative QCD predictions. These com- 
bined results are in agreement with the QCD expectations and disfavor the hypothesis 
of flavor-independent fragmentation. 

1. Introduction 

Heavy quark (Q=c,b) systems provide important laboratories for experimental tests of 
the theory of strong interactions, quantum chromodynamics (QCD). Since the large 
quark mass Mq acts as a cutoff for soft gluon radiation, some properties of these 
systems can be calculated accurately in perturbative QCD. In other cases, however, 
where QCD calculations assume massless quarks, the products of heavy hadron decays 
can complicate the comparison of data with the predictions for massless partons. It 
is therefore desirable to measure properties of both light- and heavy-quark systems as 
accurately as possible. 

In this paper we consider one of the most basic observable properties of high energy 
particle interactions, the multiplicity of charged hadrons produced in the final state. 
We consider hadronic Z° decays, which are believed to proceed via creation of a primary 
quark- ant iquark pair, Z° — > qq, which subsequently undergoes a fragmentation process 
to produce the observed jets of hadrons. If the primary event flavor q can be identified 
experimentally, one can measure the average charged multiplicity n q in events of that 
flavor, for example q — b,c, uds, where uds denotes the average over events of the types 
Z° — > mm, dd, and ss. These are not only important properties of Z° decays, but, if 
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the average decay multiplicity of the leading hadrons that contain the primary heavy 
quark or antiquark is subtracted from Uq to yield the average non-leading multiplicity, 
can also be used to test our understanding of the quark fragmentation process and its 
dependence on the quark mass. The hypothesis of flavor-independent fragmentation 
|], 0| implies that this non-leading multiplicity in e + e~ — > QQ ("heavy quark") events 
at center-of-mass (cm.) energy W should be equal to the total multiplicity in e + e~ — > 
uu,dd, and ss ("light quark") events at a lower cm. energy given by the average 
energy of the non-leading system, E nl = (1 — (xe q ))W, where (xe q ) = 2{Eq) /W is 
the mean fraction of the beam energy carried by a heavy hadron of flavor Q. 

Perturbative QCD predictions have been made H of the multiplicity difference 
between heavy- and light-quark events, Auq = Uq — n u ds- In this case the suppression 
of soft gluon radiation caused by the heavy quark mass leads to a depletion of the non- 
leading multiplicity, and results in the striking prediction that Auq is independent 
of W at the level of ±0.1 tracks. Numerical predictions of An^ = 5.5 ± 1.3 and 
An c = 1.7 ± 1.1 were also given ||. More recently, improved calculations have been 
performed [|J, confirming that the energy- dependence is expected to be very small and 
predicting Arl b =3.53±0.23 and An c =1.02±0.24 at W = M z o. 

In our previous paper || we measured n& and An& using the sample of about 10,000 
hadronic Z° decays recorded by the SLD experiment in the 1992 run. By comparing 
with similar measurements at lower cm. energies [1], |6], [7], ^| we found that An& was 
consistent with an energy-independent value, and in agreement with the prediction 
of 0. This result was subsequently confirmed by the DELPHI || and OPAL []X0f| 
Collaborations. The dominant uncertainty in our measurement resulted from lack of 
knowledge of the charged multiplicity in Z° — > cc events, n c . In this paper we present 
simultaneous measurements of rib, n c and n u d s based upon the sample of about 160,000 
hadronic Z° decays collected by SLD between 1992 and 1995, and using the SLD micro- 
vertex detector and tracking system for flavor separation. By measuring n c and n u ds 
directly we have reduced the systematic uncertainty on An b substantially, and have 
also derived An c , which allows us to compare with the QCD predictions for the charm 
system and with the only other measurement of this quantity [[l(| at the Z° resonance. 
This measurement supersedes our previous measurements of nb and Afib f|. 
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2. Apparatus and Hadronic Event Selection 



The e + e annihilation events produced at the Z° resonance by the SLAC Linear Col- 
lider (SLC) were recorded using the SLC Large Detector (SLD). A general description 



of the SLD can be found elsewhere |TTJ . The trigger and selection criteria for isolating 
hadronic Z° boson decays are described elsewhere [p~2|] . 

The analysis presented here used the charged tracks measured in the central drift 
chamber (CDC) Jl3| and in the vertex detector (VXD) [Q. A set of cuts was applied 
to the data to select well-measured tracks, which were used for multiplicity counting, 
and events well-contained within the detector acceptance. The well-measured tracks 
were required to have (i) a closest approach transverse to the beam axis within 5 cm, 
and within 10 cm along the axis from the measured interaction point; (ii) a polar 
angle 9 with respect to the beam axis within | cos 9 |< 0.80; and (iii) a momentum 
transverse to the beam axis p± > 0.15 GeV/c. Events were required to have (i) a 



minimum of seven such tracks; (ii) a thrust axis [15| direction within | cos 6^ |< 0.71; 
and (iii) a total visible energy E vis of at least 20 GeV, which was calculated from the 
selected tracks assigned the charged pion mass; 114,499 events passed these cuts. The 
background in the selected event sample was estimated to be 0.1 ± 0.1%, dominated 
by Z° — > t + t~ events. 

While the multiplicity measurement relied primarily on information from the CDC, 
the additional information from the VXD provided the more accurate impact parameter 
measurement, and D meson vertex reconstruction, used for selecting samples enriched 
in light (u,d,s) and b events, and c events, respectively. In addition to the requirements 
for well-measured tracks, "impact parameter quality" tracks were required to have (i) 
at least one VXD hit; (ii) a closest approach transverse to the beam axis within 0.3 
cm, and within 1.5 cm along the axis from the measured interaction point; (iii) at least 
40 CDC hits, with the first hit at a radius less than 39 cm; (iv) an error on the impact 
parameter transverse to the beam axis less than 250 /im; and (v) a fit quality of the 
combined CDC+VXD track x 2 /d-o.f < 5. We also removed tracks from candidate 
and A decays and 7-conversions found by kinematic reconstruction of two-track 
vertices. 

All impact parameters used in this analysis were for tracks projected into the (x—y) 
plane perpendicular to the beam axis, and were measured with respect to an average 
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primary vertex. The average primary vertex was derived from fits to ~30 sequential 
hadronic events close in time to the event under study, with a measured precision of 
&pv — (7± 2)yum [[l(J. The impact parameter 8 was derived by applying a sign to the 
distance of closest approach such that 8 is positive when the vector from the primary 
vertex to the point at which the track intersects the thrust axis makes an acute angle 
with respect to the track direction. Including the uncertainty on the average primary 
vertex the measured impact parameter uncertainty for the overall tracking system 



approaches 11 /zm for high momentum tracks, and is 76 /xm at p±Vsin9=l GeV/c Ul6| . 

3. Selection of Flavor- Tagged Samples 

We divided each event into two hemispheres separated by the plane perpendicular to the 
thrust axis. We then applied three flavor tags to each hemisphere. In order to reduce 
potential tagging bias we measured the average charged multiplicity in hemispheres 
opposite those tagged. Impact parameters of charged tracks were used to select enriched 
samples of b or light quark hemispheres, and reconstructed charmed mesons were used 
to select c quark hemispheres. 

In each hemisphere we counted the number of impact parameter quality tracks n sig 
that had an impact parameter significance of 8 norm = 8 /as >3.0. Fig. 1 shows the 
distribution of n sig upon which is superimposed a Monte Carlo simulated distribution in 
which the flavor composition is shown. For our Monte Carlo study we used the JETSET 



7.4 event generator 117] with parameter values tuned to hadronic e + e annihilation data 



IS], combined with a simulation of 5-decays tuned to T(4S) data and a simulation 



of the SLD. A more detailed discussion of flavor tagging using impact parameters can 
be found in |16| . The Monte Carlo simulation reproduces the data well and shows that 
most light quark hemispheres have n S i g =0 and that the n sig >3 region is dominated 
by b quark hemispheres. Hemispheres were tagged as light or b quark by requiring 
n S i g = or n sig > 3, respectively. Table 1 shows the number of light and b quark 
tagged hemispheres and their flavor compositions estimated from the simulation. 

From Fig. 1 it is clear that an impact parameter tag does not provide a high- 
purity sample of c quark hemispheres. For this purpose we required at least one 
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uds-tag 


c-tag 


6-tag 


jfc hemispheres 


154,151 


976 


9,480 


composition 


uds 


0.752±0.001 ±0.004 


0.074±0.002 ±0.014 


0.014±0.001 ±0.001 


c 


0.158±0.001 ±0.006 


0.640±0.008 ±0.025 


0.048±0.001 ±0.005 


b 


0.089±0.001 ±0.004 


0.286±0.005 ±0.022 


0.938±0.001 ±0.006 



Table 1. Numbers of hemispheres and fractional compositions of uds, c and b quarks in the tagged 
hemispheres. The first quoted errors represent the errors due to the limited size of the Monte Carlo 
sample and the second are due to the uncertainties from the modelling of heavy hadron production 
and decay. 



prompt D* + or D + mesonQ reconstructed in a hemisphere. This tag is similar to that 
described in |L9" |. The D* + mesons were identified using the decay D* + — > tt^D , where 
7r+ is a low-momentum pion and the D° decays via D° — > K~n + ("three-prong"), 
D° -> K-it + it ("satellite"), or D° -> R-tt+tt+tt- ("five-prong") modes. The D+ 
mesons were indentified using the decay mode D + — > K~tt + tt + . D meson candidates 
were formed from all combinations of well- measured tracks with at least one VXD hit. 
D° candidates were formed by combining two (for the three-prong and satellite modes) 
or four (for the five-prong mode) charged tracks with zero net charge, and by assigning 
the K~ mass to one of the particles and tt + mass to the others. 

For D* + candidates, we first required a candidate D° in the mass range 1.765 
GeV/c 2 < M™ nd - <1.965 GeV/c 2 (three-prong), 1.815 GeV/c 2 < Mg* nd <1.915 
GeV/c 2 (five-prong), or 1.500 GeV/c 2 < Mffi d - < 1.700 GeV/c 2 (satellite). D* + candi- 
dates were then required to pass either a set of kinematic cuts or a set of decay length 
cuts to suppress combinatorial backgrounds and backgrounds from B — ► D* + decays. 
The kinematic cuts are: (i) \ cos6kd°\ < 0.9 (three-prong and satellite modes) and 
| cos9jcd°\ < 0.8 (five-prong mode), where 6kd° is the angle between the D° direction 
in the laboratory frame and the K direction in the D° rest frame, (ii) p n + >1 GeV/c, 
and (iii) xe d „ + >0.4 for the three-prong and satellite modes and XE Dt+ >0.5 for the 
five-prong mode, where xe d * + = 2Ed*+/W and Ee>*+ is the D* + energy. For the decay 
length analysis we performed a fit of the D° tracks to a common vertex and calculated 
the decay length, L°, between the primary vertex and this D° decay vertex, and its 



1 In this paper charge- conjugate cases are always implied. 
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error, a L o. The decay length cuts are: (i) a \ 2 probability>l% for the vertex fit to the 
D° tracks, (ii) a decay length significance L° /a L o >2.5, (iii) the two-dimensional im- 
pact parameter of the D° momentum vector to the interaction point <20/zm, and (iv) 
x E Dt+ >0.2 for the three-prong and satellite modes and xe d * + >0.4 for the five-prong 
mode. 

For all D* + candidates we required the proper decay time of the D°, r propei . = 
L° / (3^1 — /3' 2 , where (3 = pi>o/E D o and pr>o and E D o are the reconstructed momentum 
and energy, respectively, of the candidate D° meson, to be in the range 0< r proper <lps. 
Figs. 2(a), (b) and (c) show the distribution of AM, where AM = M^ t n + d - M^ d , after 
the above cuts for the three D° decay modes, upon which is superimposed the Monte 
Carlo simulated distribution in which the flavor composition is shownP]. A hemisphere 
was tagged as c if it contained a D* + candidate with AM <0.15 GeV/c 2 . 

D + — > K~7t + tt + candidates were formed by combining two tracks of the same sign 
with one track of the opposite sign, where all three tracks were required to have momen- 
tum p >1 GeV/c. The two like-sign tracks were assigned n + masses, the opposite-sign 
track was assigned the K~ mass, and all three tracks were fitted to a common ver- 
tex. A series of cuts was applied to reject random combinatoric, D* + , and B-decay 
backgrounds. We required: (i) xe d+ > 0.4, (ii) cos9kd+ > —0.8, where 6kd+ is 
the angle between the directions of the D + in the laboratory frame and the K~ in 
the D + rest frame, (iii) the mass differences M(K~-k + 'k + ) — M(K~tt + ) for each of 
the two pions to be greater than 0.16 GeV/c 2 , (iv) the normalized D + decay length 
L + /cr^+ >3.0, and (v) the projection of the angle between the D + momentum vec- 
tor and the vertex flight direction to be less than 5 mrad in the plane perpendicular 
to the beam axis and less than 20 mrad in the plane containing the beam axis. A 
hemisphere was c-tagged if it contained a D + — > K~h + 'k + candidate in the mass range 
1.800 GeV/c 2 < M(K-n + 7i + ) <1.940 GeV/c 2 . Fig. 2(d) shows the mass M(K-n + 7i + ) 
distribution of the data upon which is superimposed the Monte Carlo simulated dis- 
tribution in which the flavor composition is shown. 

The union of the three samples of D* + candidates and the sample of D + candi- 

2 In the Monte Carlo simulation the production cross section and branching fractions, and nor- 
malization of the AM distributions in the region AM > 0.15 GeV/c 2 , for the D° — » K~n + and 
D° K~ir + ir° modes were adjusted to match the data in Fig. 2, as described in Ref. [19]. The 
adjustment was small and included in the systematic errors. 
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dates was used to tag c quark hemispheres. The flavor composition of these tagged 
hemispheres is shown in Table 1. Approximately 400 of the c-tagged hemispheres were 
also tagged as either b or uds hemispheres. Monte Carlo studies indicated that these 
were mostly true c-hemispheres. The exclusion of these hemispheres from the b- and 
uds-tagged samples was found to have negligible effect on the final results. 



4. Measurement of Charged Multiplicities 

Well-measured charged tracks defined in Section 2 were counted in the hemispheres op- 
posite those tagged. The measured average hemisphere multiplicities m; (i = uds, c, b) 
were m u d s — 8.94 ± 0.01, m c = 9.15 ± 0.12 and m& = 9.99 ± 0.04 (statistical errors 
only). 

The rrii are related to the true average multiplicities nj (j = uds, c, b) of uds, c and 
b quark events by: 

2xm, = Pi Ms C iMs n uds + P ijC (C*nf + Cf c n n J) + + C f W) (1) 



where: Pij is the fraction of hemispheres of quark type j in the i-tagged hemisphere 
Tj = nj k + nj 1 (j 7^ uds), and Itf 



sample; n,- = n^ k + rf^ 1 (j ^ uds), and n^ k is the true average multiplicity originating 



from the decay products of j-hadrons and rij 1 is that originating from the non-leading 
particles; C ijU d s is the ratio of the average number of measured charged tracks in light 
quark hemispheres opposite z-tagged hemispheres, to the average number of charged 
tracks in true light quark hemispheres; Cf k (j ^ uds) is the ratio of the average num- 
ber of measured charged tracks originating from the decay products of j-hadrons0 in 
hemispheres opposite i-tagged hemispheres, to the average number of tracks originating 
from the decay products of j-hadrons; C"j (j ^ uds) is the ratio of the average number 
of measured charged tracks originating from the non-leading particles in true j-quark 
hemispheres opposite those tagged as i-quark hemispheres, to the average number of 
tracks originating from non-leading particles in true j-quark hemispheres. The con- 
stants P are shown in Table 1. The constants C were also calculated from our Monte 
Carlo simulation and are shown in Table 2; they account for the effects of detector ac- 
ceptance and inefficiencies, for tracks from beam-related backgrounds and interactions 



3 We include the products of both strongly and weakly decaying heavy hadrons. 
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3 


uds 


c 


b 


i 


dk 


nl 


dk 


nl 


uds 


0.875±0.001±0.001 


0.798±0.002±^ 


0.885±0.002lg;g^ 


0.820±0.002±S 


0.887±0. 003±g;g|? 


c 


0.803±0-019±g:gg| 


0.831±0.011±0.004 


0.864±0. 0091^17 


0.854±0.013+^;^4 


0.849±0.015±0.026 


b 


0.875±0.015t°;°°^ 


0.816±0.013±0.003 


0.887±0.010i°;ol8 


0.854±0.003^°;°2i 


0.893±0.004±0.025 



Table 2. The constants C calculated from the Monte Carlo simulation. The first quoted errors are 
statistical and arise from the finite size of the Monte Carlo sample. The second are due to the 
uncertainties from C and B hadron production and decay. 



in the detector material, and for biases introduced by the event and tagged-sample 
selection criteria. We included in the generated multiplicity any prompt charged track 
with mean lifetime greater than 3 x 10~ 10 s, or any charged decay product with mean 
lifetime greater than 3 x 10 _10 s of a particle with mean lifetime less than 3 x 10~ 10 s. 
We fixed nf — 5.20 and =11.10, using the measured values from [^, ^T], 22] with 



the addition of 0.20 and 0.22 tracks, respectively, estimated from the Monte Carlo 
simulation, to account for the effects of higher mass states of heavy hadrons produced 
in Z° decays. 

We then solved eqns. (1) to obtain the average charged multiplicities per event, 
n uds = 20.21 ± 0.10, n c = 21.28 ± 0.46 and n b = 23.14 ± 0.10 (statistical errors only). 
The multiplicity differences between c and light quark events, and b and light quark 
events are, respectively 

An c = 1.07 ±0.47 (stat.) 
An b = 2.93 ±0.14 (stat.). 

5. Systematic Errors 

Experimental systematic errors arise from uncertainties in modelling the acceptance, 
efficiency and resolution of the detector. Systematic uncertainties also arise from er- 
rors on the experimental measurements that function as the input parameters to the 
modelling of the underlying physics processes, such as errors on the modelling of b and 
c fragmentation and decays of B and C hadrons. 

The effect of uncertainty in the tracking efficiency was estimated to cause a common 
±0.9% variation of the constants C. The effect of uncertainty in the corrections for the 
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Source of Uncertainty 


rinds n c n b An c An b 


Tracking efficiency 
7 conversion & fake tracks 
Monte Carlo statistics 


±0.182 ±0.194 ±0.205 ±0.012 ±0.023 
±0.101 ±0.108 ±0.114 ±0.007 ±0.013 
±0.046 ±0.212 ±0.045 ±0.217 ±0.064 


Total 


±0.213 ±0.307 ±0.239 ±0.217 ±0.069 



Table 3. Systematic errors due to detector modelling. 



residual 7 conversions and fake tracks was estimated to cause a common ±0.5% vari- 
ation of the constants C . Statistical effects from the limited Monte Carlo sample size 
were also considered. These errors, summarized in Table 3, were added in quadrature 
to obtain a total systematic error due to detector modelling. Note that the uncertain- 
ties in total track reconstruction efficiency are the dominant source of systematic error 
for n u ds and n b , but are small for the differences An c and An b . In the case of n c , An c 
and An b the statistical error from the limited Monte Carlo sample size is dominant. 

We performed several consistency checks on our results. We checked that our 
Monte Carlo simulation showed good agreement with the data for track p± and cos 9 
distributions in the hemispheres opposite those tagged. We then varied the thrust 
axis containment cut within 0.5< |cos#t| <0.8. To check for possible bias from our 
hemisphere tags the cut on the track significance S norm was varied from 2.0 to 4.0 
for the light and b quark hemisphere tags, and D* + and D + mesons were considered 
separately as a c quark hemisphere tag. We also removed hemipheres tagged as both 
c and uds or h. Finally, we performed our analysis separately in 2- and > 3-jet event 
samples selected using the Durham algorithm |23| with y cut =0.003, to check for any 
possible bias in multi-jet events. In each case all the re-evaluated Hi were found to be 
consistent with our central values of Hi within the statistical errors. 

In order to estimate the systematic errors due to uncertainties in modelling heavy 
hadron production and decay we used an event re-weighting scheme to vary the mul- 
tiplicity distributions in the Monte Carlo simulation and to obtain modified values of 
the constants C and P. The effect of uncertainty in heavy flavor fragmentation was 



estimated by varying the e parameter of the Peterson fragmentation function [24] used 
as input to generate the Monte Carlo sample, corresponding to 5 < xe >=±0.012 and 
±0.011 for c and b quarks respectively, corresponding to the average errors in measure- 
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Source of Uncertainty 


Variation 




n c 


n b 


An c 


An b 




b fragmentation 


(x Eb )=0.700±0.011 


±0.001 


+0.002 
-0 


+0.288 
-0.281 


+0.004 

-0 


+0.289 
-0.279 




B meson lifetime 


t 6 =1.55±0.1 ps 


±0.001 


+0.027 
-0.028 


+0.010 
-0.007 


+0.027 
-0.026 


+0.012 
-0.007 




B baryon lifetime 


r 6 =1.10±0.3 ps 


n 
+u 

-0.008 


-0.036 


i n nn8 

TU.UUo 

-0.001 


-0 


-\- u . U 1 z 

-0 




B baryon prod, rate 


/a„ = 9% ± 3% 


+0.004 
-0.001 


+0 

-0.001 


+0.021 
-0.020 


+0.001 
-0.003 


+0.019 
-0.018 




Rb (b fraction) 


0.221±0.003 


±0.001 


+0 007 
-0.006 


+0 041 
-0^040 


+0 007 
-0^008 


+0 040 
-0.039 




B -> D+ +X fraction 


0.17±0.06 


+0 

-0.007 


+0.054 
-0 


+0 

-0.036 


+0.053 

-0 


+0 

-0.024 




c fragmentation 


(£,e c )=0.494±0.012 


+0.008 
-0.010 


+0.236 
-0.155 


+0.004 
-0.002 


+0.244 
-0.151 


+0.015 
-0.008 




R c (c fraction) 


0.l7l±0.020 


+0.026 
-0.027 


+0.081 
-0.099 


+0.006 
-0.007 


+0.107 
-0.126 


±0.033 




cc — > D + + X fraction 


0.20±0.04 


+0.004 
-0.003 


+0.035 
-0.039 


±0.006 


+0.039 
-0.042 


+0.010 
-0.009 




Tprdk 
ll c 


5.20±0.26 


±0.003 


+0.010 
-0.009 


+0.001 

-0 


+0.005 
-0.006 


±0.003 




n b 


I1.10±0.36 


±0.003 


+0.009 
-0.008 


±0.016 


±0.012 


±0.013 


D°- 


■+ K~ir + , D° — > A'~7r + 7r° production 


-20% 


-0.013 


+0.062 


-0.003 


+0.075 


+0.010 


Total 




+0.028 
-0.034 


+0.269 
-0.194 


+0.293 
-0.287 


+0.289 
-0.203 


+0.296 
-0.286 



Table 4. Systematic uncertainties due to heavy hadron modelling. 



ments of these quantities p5| . The average B hadron lifetime was varied by ±0.1 ps for 
B mesons and ±0.3 ps for B baryons p6 |. The effect of varying the B baryon production 
rate in b events by ±3% |16| was also examined. Absolute variations of ±6% and ±4% 
were applied to the B — > D + branching ratio and c — > D + branching ratio, respectively 
16| . The effect of the present experimental uncertainties in the branching fractions, 



R c = T(Z° -> cc)/T(Z° -> qq) and R b = T(Z° -> bb)/T(Z° -> qq), of 5R C =±0.020 and 
<5_Rf,=±0.003 respectively |^6| were also included. The decay multiplicities of C and B 
hadrons were varied by ±0.26 and ±0.36 charged tracks, respectively |20], |21], [Z2|. For 
the D* + analysis we also accounted for the adjustment of the production cross section 
and branching fractions for the D° — > K~ir + and D° — > K~7t + tt° modes in the Monte 
Carlo by assigning the full shift of the Monte Carlo simulated distribution as a sys- 
tematic error. These uncertainties, summarized in Table 4, were added in quadrature 
to obtain total systematic uncertainties due to C and B hadron modelling. For An c 
(An;,) the dominant contributions were from the uncertainties in c (b) fragmentation 
and R c (Rb). 
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6. Summary and Conclusions 



Combining systematic uncertainties in quadrature we obtain: 

n uds = 20.21 ± 0.10 (stat.) ± 0.22 (syst.) 
n c = 21.28 ± 0.46 (stat.) + (syst.) 
n b = 23.14 ± 0.10 (stat.) + (syst.). 

Subtracting nf fc =5.20 and n^ fc =11.10 from our measured n c and n b respectively, we ob- 
tained the average non-leading multiplicities rtf = 16.08 ± 0.46(stat.) la^syst.) an d 
rtfi 1 = 12.04 ± 0.10(stat.) l^l^syst.). The hypothesis of flavor-independent fragmenta- 
tion implies that TIq(W) = n uds ([l — (xe q (W))]W). Fig. 3(a) shows our measurement 
of n u ds plotted at W = Mz, and our measurements of nf and n b l plotted at the ap- 
propriately reduced non-leading energy [1 — (x Eq (W))]W . Previous measurements of 
these quantities |], §, |, |, 0, ^ IH 

are also shown. The curve is a fit to the energy 
dependence of the n u d s measurements shown and those at 5 < W < 92 GeV [p7| . Fig. 
3(b) shows the differences between the non-leading data points in Fig. 3(a) and the 
curve. A linear fit to these differences (Fig. 3(b)) yields a slope of s = 1.14 ± 0.32 
tracks/ln(GeV). This differs from the expectation for identical energy dependence, 
s = 0, by 3.6 standard deviations, indicating that the hypothesis of flavor-independent 
fragmentation is disfavored at this level. 

Combining systematic uncertainties in quadrature we obtain: 

An c = 1.07 ± 0.47 (stat.)+ (syst.) 
An b = 2.93 ± 0.14 (stat.)+ (syst.). 

Fig. 4 shows our measurements of An c and An b together with those from other ex- 
periments |I1 |5], |7j, ^, [TTJ, E7J, at the respective cm. energies. The new result for 
An b is consistent with our previous measurement || and with the measurements from 



LEP |9|, [Uj and Mark-II pTl , and that for An c is consistent with the OPAL measure- 



ment [|Kj]. Linear fits to the An c and An b data as a function of ln(VK) yield slopes 
of s=— 1.33±1.04 and s=— 1.43±0.82 tracks/ln(GeV), respectively. These slopes are 
consistent with the perturbative QCD prediction of energy independence [3], s=0, at 
the level of 1.3 and 1.7 standard deviations, respectively. 

Comparing our measurements of An c and An b with the predictions of Refs. [3,4] 
(Fig. 4) we found that both were in good agreement with the predictions of Ref. [4], 
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while the former was in good agreement with the prediction of Ref . [3] , and the latter 
within 1.7a of this prediction. 

As a result of the accurate measurements of An c and An h at W = M z o , constraints 
on the energy dependence of these quantities are now limited by the uncertainties in 
the lower energy measurements. In order to improve the constraints on the validity of 
perturbative QCD calculations at the scales M b or M c , it is necessary to improve the 
accuracy of the measurements of An h and An c , respectively, at lower energies, and/or 
extend the ln(W) lever-arm of such measurements. It would thus be desirable to have 
measurements of An c from the continuum below the T(4S), and for both An c and An^ 
to be measured at LEP-II and e + e~ colliders at even higher energies. 
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Figure captions 



Figure 1. The distribution of the number of tracks per hemisphere n S i g that miss 
the interaction point by more than 3a in the x-y plane. The points represent the data 
distribution and the solid histogram represents the Monte Carlo simulated distribution. 
The flavor composition of the Monte Carlo distribution is shown. 

Figure 2. The distributions of AM for a) D° -> K~ir + , b) D° -> R-n+n and c) 
D° -> K-7T+7T+7T-; d) M(K-iT + 7T + ) distribution for D + -> K~7r + 7r + (see text). The 
points represent the data distributions and the solid histograms represent the Monte 
Carlo simulated distributions. The flavor composition of the Monte Carlo distributions 
is shown. 

Figure 3. a) Our measurements of n u ds plotted at W = M^o and the non-leading 
multiplicities n^ 1 and n£ l plotted at the appropriately reduced non-leading energy [1 — 
(xe q {W))]W . Previous measurements of these quantities |], || |7|, || ||, |H], 
also shown. The solid curve is a fit |27j to n u d s measured in the range 5 < W < 92 GeV. 
The error on this curve (dotted lines) is dominated by the uncertainty on the removal 
of the heavy quark contribution to each measured total charged multiplicity, b) The 
differences (points) between the non-leading data points in a) and the solid curve. A 
linear fit to these differences is shown by the dashed line. For clarity the different data 
points at the same energy are displayed with small relative displacements in W. 

Figure 4. Multiplicities differences a) An c and b) An?, as functions of cm. energy. 
The predictions of Ref. [3] are shown as the solid lines and those of Ref. [4] are 
shown as the dashed lines. For clarity the different data points at the same energy are 
displayed with small relative displacements in W. 
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